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Abstract 



Introduction shape measurement: One is the circularization of the im- 



o , 

^\i I We examine the anisotropic point spread function (PSF) of Suprime-Cam data utihzing the dense star 

field data. We decompose the PSF eUipticities into three components, the optical aberration, atmospheric 
^ i' turbulence and chip-misalignment in an empirical manner, and evaluate the amplitude of each component. 

Then, we test a standard method for correcting the PSF eUipticities used in weak lensing analysis against 
mock simulation. We find that, for long-exposure data, the optical aberration has the largest contribution 
to the PSF eUipticities, which can be modeled well by a simple analytic function based on the lowest-order 
aberration theory. Statistical properties of PSF eUipticities originated from the atmospheric turbulence are 
investigated by using the numerical simulation, and it is found that simulation results are in a reasonable 
agreement with the observed data. It follows from those findings that the spatial variation of PSF eUiptic- 
ities consists of two components; one is a smooth and parameterizable component arisen from the optical 
PSF, and the other is a non-smooth and stochastic component originated from the atmospheric PSF. The 
O,' former can be well corrected by the standard correction method with polynomial fitting function. However, 

I , for the later, its correction suffers from the common limitation arisen from sparse sampling of PSFs due 

2 ' to a limited number of stars. We also examine effects of the residual PSF anisotropics on Suprime-Cam 

*— > I cosmic shear data (5.6-degree^ of i'-band data). We find that the shape and amplitude of B-mode shear 

jrt ' variance are broadly consistent with those of the residual PSF eUipticities measured from the dense star 

field data. This indicates that most of the sources of residual systematic are understood, which is an 
important step for cosmic shear statistics to be a practical tool of the precision cosmology. 
►^ Key words: cosmology: observations — dark matter — large-scale structure of universe 

—L. ' ages which systematically lower the lensing shear. The 
• , Weak gravitational lensing has now became a unique other is a coherent deformation of the images caused by 
^i ' and practical tool to probe dark matter distribution ir- anisotropy of the PSF which may mimic a weak lensing 
Z^ i respective of its relation to luminous objects, thanks to shear. A lot of effort has been made to develop tech- 
great progress in technique for weak lensing analysis as niques and software to correct the PSF effects, and the 
well as in wide field imaging instrument (see for reviews weak lensing community conducted blind tests using mock 
Mellier 1999; Bartelmann & Schneider 2001; Refregier data to evaluate their performance (Hcymans et al. 2006a; 
2003; Hoekstra & Jain 2008; Munshi et al. 2008). Massey et al. 2007; Bridle et al. 2010; Kitching et al. 
;h ■ Measurement of statistics of weak lensing shear, also 2012a, b; and references therein). Most of the softwares 
^.' known as the cosmic shear statistics, has been recog- are designed to correct the PSF regardless of its origins, 
nized as one of most powerful probes of the cosmologi- Major sources of the PSF, recognized so far except for the 
cal parameters and is employed as the primary science diffraction, include the atmospheric turbulence, the opti- 
application by many on-going/future wide field survey cal aberration, the misalignment of CCD chips on a focal 
projects (e.g., the Panoramic Survey Telescope & Rapid plane and the pixelization. Alternative and complemcn- 
Response System (Pan-STARRS); Dark Energy Survey tary approach to this issue is to investigate properties of 
(DES); Hyper SuprimeCam survey; the Kilo-Degree the PSF for a specific instrument with focusing on a spe- 
Survey (KIDS); the Large Synoptic Survey Telescope cific cause(s), which helps to optimize the PSF correction 
(LSST); and Euclid mission.) It should be, however, scheme and to understand a level of residual systemat- 
noted that the full potential of the cosmic shear statis- ics (Hoekstra 2004; Jarvis & Jain 2004; Wittman 2005; 
tics to place a constraint on the cosmological parameters Jarvis, Schechter & Jain 2008; Jee et al. 2007; 2011; 
is archived only if systematic errors are reduced down to Rhodes et al. 2007; Jee & Tyson 2011; Heymans et al 
a required level (Huterer ct al. 2006). 2012; Chang ct al. 2003). 

In the process of weak lensing analysis, the correction The latter approach is exactly what we take in the 

for the point spread function (PSF) is of fundamental present paper. The purpose of this study is twofold: 

importance, because it has two serious effects on galaxy First, we investigate properties of anisotropic PSFs of 
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Subaru prime focus Camera (Suprime-Cam, Miyazaki et 
al. 2002). We pay a special attention to evaluate the am- 
plitude and spatial correlation of PSF anisotropics caused 
by different sources. To do this, we use a set of dense star 
field images from which we can sample the PSF densely 
and investigate the spatial variation of the PSF in detail. 
Second, using mock simulation data of PSFs, we examine 
how well the PSF anisotropy can be corrected with the 
standard PSF correction scheme. 

The motivation of this study is the existence of non-zero 
B-mode in the cosmic shear correlation functions mea- 
sured from Suprime-Cam data. The measurement was 
made by Hamana et al. (2003), they analyzed a 2 degree^ 
data and detected statistically significant B-mode corre- 
lations (see §5 for an updated but similar result). Since 
the B-mode shear is not produced by gravitational lens- 
ing in the standard gravity theory at the lowest order of 
the perturbative treatment^, the existence indicates any 
non-lcnsing process(es) taking place. The origin(s) of the 
B-mode was unclear; it could be simply due to remain- 
ing systematics in the data reduction and analysis. Also 
it may arise from an intrinsic alignment of galaxy shapes 
(e.g.. Croft & Metzler 2000; Catelan, Kamionkowski & 
Blandford 2001; Crittenden et al. 2002; Jing 2002; Hirata 
& Seljak 2004; Heymans et al. 2006b). More importantly, 
the B-mode may be produced by a non-standard grav- 
ity. Thus the B-mode shear, in principle, provides us with 
valuable information on the gravity theory. It is thus im- 
portant to first understand the level of residual systematic 
in the data analysis. 

The outline of this paper is as follows; in Section 2, 
data reduction and PSF measurement of dense star field 
data is described. The amplitude and spatial correlation 
of PSF anisotropics caused by different sources are evalu- 
ated separately in Section 3. The PSF correction scheme 
adopted in our cosmic shear analysis is tested against the 
mock simulation data in Section 4. In section 5, we assess 
the impact of the residual PSF anisotropy on the cosmic 
shear analysis. Finally, summary and discussion are given 
in Section 6. In Appendix, we summarize the proper- 
ties of PSF ellipticitics caused by the third-order optical 
aberrations (Appendix 1), and a numerical study of PSF 
ellipticitics caused by atmospheric turbulence (Appendix 
2). 

2. Data reduction and PSF measurement of 
dense star field images 

We use i'-band data of a dense stellar field taken with 
the Suprime-Cam on 2002/9/30, 2003/6/30 and 2003/7/1. 
The field is located at the Galactic coordinate of {l,b) = 
(38°,— 3°). We collected 289 shots of 60 seconds exposure 
and 70 shots of 30 seconds exposure from the data archive 
SMOKA^. The seeing FWHMs of those data range from 
0.46"to 0.88"with the median of 0.63". 



Each CCD data was reduced using the SDFred'^ soft- 
ware (Yagi et al. 2002; Ouchi et al. 2004). Note that 
we conservatively use the data only within 15 arcmin ra- 
dius from the field center of Suprime-Cam, because at 
the outside of that, the point spread function (PSF) be- 
comes elongated significantly which may make the cor- 
rection for the PSF inaccurate. Then mosaicking of 10 
CCDs was performed with SCAMP (Bertin 2006) and 
SWarp"^ (Bertin et al. 2002). In addition to the individ- 
ual exposure images, we also generate stacked images us- 
ing SWarp. Note that the images obtained during a same 
night arc taken with the same pointing (thus no dithering 
was made) . Since differences in pointings of images taken 
in different nights are very small (less than the gaps of 
CCDs), no data from different CCDs is added. 

Object detections were performed with SExtractor 
(Bertin & Arnouts 1996) and hfindpeaks of IMCAT soft- 
ware (Kaiser, Squires & Broadhurst 1995), and two cat- 
alogs were merged by matching positions with a toler- 
ance of 1 arcsec. From the image of stars, we measure 
an anisotropy of the PSF. We use the following ellipticity 
estimator (exactly speaking the polarization, see Kaiser 
et al. 1995) to characterize the anisotropy of PSFs: 
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where f{6) is the surface brightness of an object and 
Wcid) is the Gaussian window function. We used get- 
shapes of IMCAT for the actual computation of these 
quantities. 

From catalogs of detected objects (almost all of them 
are stars) , we generate two kinds of sub-catalogs by mim- 
icking weak lensing analysis: One is "star-role" catalog 
which contains 700 (the number density of 1 arcmin^^) of 
the highest SN non-saturated stars. The other is "galaxy- 
role" catalog which contains 30,000 (42 arcmin"^) of the 
next highest SN non-saturated stars. The minimum flux 
SNs for star-role and galaxy-role catalogs are sufficiently 
high; SNmin > 330 and > 40, respectively. 

3. Analysis of PSF ellipticities 

An anisotropy of PSF is resulted from several sources 
including the atmospheric turbulence, the optical aberra- 
tion, the misalignment of CCD chips on a focal plane, the 
pixelization, and tracking error of the telescope. Our at- 
tempt in this section is to examine the PSF on a source- 
by-source basis using physically motivated models. We 
evaluate the statistical properties of PSF anisotropics 
using a two-point correlation function related estimator 



The B-mode shear can arise from higher order terms but their 
expected amphtude is much smaller than the observed signal 
(Schneider, van Waerbeke & Mellicr 2002; Hilbcrt et al. 2009). 
http : //smoka. nao .ac.jp/ 



^ In the process of the correction for both the field distortion 
and differential atmospheric dispersion, the bi-cubic resampling 
scheme was implemented to suppress the aliasing effect (Hamana 
& Miyazaki 2008). 

* SWarp was modified so that it can treat the bad pixel flag from 
the SDFred software properly. 
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called the E- and B-niode aperture mass variances de- 
fined by (Schneider et al. 1998; we follow the notation by 
Schneider, van Wacrbcke & Mellier 2002) 



(AOW=2 



and 



{Mim=- -^ 



20 



e+(0)r+^ -XWT-^ 



,(4) 



where T+{(j)) and T_{(j)) 



are defined in Schneider et al. 
(2002), and ^+((^) and ^-((/)) are the two-point correlation 
functions of the ellipticities computed by 



^+{d) = Y^ (e,tej-t + e,xejx) /A^p(6') 
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In the above expressions, the summation is taken over 
all the pairs of objects with a distance within a width of 
a bin considered - Abin/2 < (f) < 9 + Abin/2, Np{9) is 
the number of pairs in the bin, and et and Cx are the 
tangential and 45° rotated ellipticity components in the 
frame defined by the line connecting a pair of objects. 

3.1. Optical aberrations 

Although the PSF ellipticities arise from several causes 
and vary temporally, there is a visibly identifiable com- 
ponent that is most likely originated from optical aberra- 
tions. In top-left panel of Fig. 1 and 2, we show two exam- 
ples of whisker plots of PSF ellipticities where characteris- 
tic features of optical aberration arc clearly appeared, and 
we also show one case in Fig. 3 where optical aberration 
is almost invisible. We fit the spatial variation of PSFs 
with the following fitting function of the optical PSF ellip- 
ticities which is motivated from the lowest-order optical 
aberration theory (see Appendix 1 for details), 
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cos2f/' 
sin 2?/; 



cos(V' — 0) 
sin('0 — 9) 



(7) 



where rig = {1,2,3,4} and ria — {1,2}. We determine 
the model parameters s„^ , a„^ and 9 by the standard 
least square method with observed PSF ellipticities. This 
model consists of the axis-symmetric, asymmetric and 
constant components, which are displayed in the bottom- 
left panel, bottom-right panel, and top-right corner of top- 
panels, respectively. Note that in the top-left panel, the 
PSF ellipticities measured from "star-role" stars but af- 
ter a constant component being subtracted are shown. It 
is clearly seen from Fig. 1 and 2 that the simple three- 
component model can well reproduce the observed PSFs. 
It may be therefore said that the optical aberration is one 
major origin of the PSF anisotropy. Note that a cause(s) 
of the constant component is uncertain. It arises not only 
from the optical aberration (the misalignment coma, see 
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Fig. 1. Ellipticity maps, so-called whisker plot, showing an 
illustrative examples of ellipticity pattern in the Suprime-Cam 
measured from a 60-sccond exposure data. The top-left panel 
shows the observed PSF ellipticites measured from high-SN 
star images ("star- role" stars), but after the mean compo- 
nent which is displayed in the top-right corner being sub- 
tracted. We fit the spatial variation of PSFs with the opti- 
cal aberration model given by eq. (7). The model consists 
of the symmetric, asymmetric and constant components. A 
sum of the first two components is shown in the top-left 
panel, in which the constant component is plotted in the 
top-left corner. The bottom-left/right panel shows the sym- 
metric/asymmetric component individually. 
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Fig. 2. Same as Fig. 1 but for another 60-second exposure. 
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Fig. 3. Same as Fig. 1 but for another 60-second exposure. 
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Fig. 4. Scatter plot showing the amplitude of the symmetric 
and asymmetric components of the optical aberration moti- 
vated PSF model for each exposure. Crosses are for 30-sec- 
ond exposures, and small filled squares are for the 60-second 
exposures. Three examples shown in Fig. 1, 2 and 3 are 
marked with the open circle, open square and open triangle, 
respectively. Dotted boxes define groups used for statistical 
analyses. 

Appendix 1) but also from, e.g., the tracking error. Since 
the constant ellipticity component does not contribute to 
the aperture mass variances, we do not consider it in the 
following discussion. 

Having found that the optical aberration is one ma- 
jor source of the PSF cUipticitics and can be well mod- 
eled with one constant and two spatially variable com- 
ponents, we introduce estimators that quantify the am- 
plitude of the symmetric and asymmetric components: 
as = sign{ely^y/'^ and a^ = {elsym)^^"^, where the RMS 
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Fig. 5. The aperture mass variance of PSF cUipticitics as a 
function aperture radius. Top and and panels are for E- and 
B-mode, respectively. The left panels show ones measured 
from observed data, whereas the right panels are from the op- 
tical aberration model. DifEercnt lines are for different groups 
defined in Fig 4, and the plotted aperture mass variances are 
average values taken over all the exposures in each group. 

scatter (e^)^/^ is computed from the model ellipticities of 
700 star-role PSFs for each exposure, and "sign" in as is 
+ or — for the tangentially (the mean tangential ellipticity 
has the plus sign) or radially elongated cases, respectively. 
The results for all the exposures arc plotted in Fig. 4. It is 
observed from the figure that there is a "V-shaped" trend 
between strengths of the symmetric component and asym- 
metric component. In the lowest-order aberration theory, 
the symmetric aberration is caused by a shift in the focal 
plane position from an ideal position, whereas the asym- 
metric one is caused by the decenter and/or tilt between 
the axes of the focal plane and the corrector (lenses) . The 
result indicates that those two deviations happen in a mu- 
tually related manner. 

Next we examine the statistical properties of the PSF 
ellipticities using the aperture mass variance, eqs. (3) and 
(4) with especially focusing on contribution from the op- 
tical aberration. To do so, we classify data using as and 
aA into 9 groups as shown in Fig. 4. We compute av- 
erage aperture mass variances for each group and show 
results in Fig. 5. Left panels of the figure show the results 
from the observed data, from which it is found that the 
E-mode variance has a larger amplitude on larger scales, 
whereas in most cases B-mode variance is smaller than the 
E-mode especially on larger scales. If we assume a power- 
law spectrum for ellipticity power spectrum^ Pe{l) oc P, 



The aperture mass variance relates to the ellipticity power spec- 
trum as (Schneider et al. 1998) {M^){e) x J dl lPe(l)I^{W) 
where I(x) is the aperture function. Thus {M'^)(9) oc 0—"—^ for 
a power-law power spectrum of Pe(0 oc I". 
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Fig. 6. EUipticity magnitude for each component is shown: 
upper and lower panels are for ei and £2, respectively. This 
is taken form the same data as Fig. 1. The left panels show 
the observed ellipticity, whereas the central panels show the 
aberration model. The color scale of the central panels are 
same as that of the left panels. The right panels show the 
residual ellipticity (cobs ~ ^opt— model)- Note that the color 
scale of the right panels is shown in just bottom of each plot, 
and are different from the other panels. 

the measured E-mode variance indicates that n ~ —3. 

Right panels of Fig. 5 show the aperture mass variances 
computed from PSF eUipticities of the optical aberration 
model. Comparing observed variance in the left panels 
with ones in the left panels, we can estimate contribu- 
tion from the optical aberration to PSF ellipticity. We 
first look at the E-mode. It follows from the comparison 
between the left and right panel of Fig. 5 that the op- 
tical aberration is the dominant component of the PSF 
eUipticities at least for large aberration cases such like 
Gl-3, G8 and G9 as expected from the visual impres- 
sion (Figs. 1 and 2). For small aberration cases (G4-G7), 
one may see the aperture mass variances of the model are 
slightly smaller than the observed ones. In the following 
subsections, we argue that the residual variances mostly 
arise from atmospheric turbulence and chip misalignment. 
Let's turn to the B-mode. It follows from the comparison 
between the left and right panel that the variances mea- 
sured from the model are smaller than the observed ones 
except for Gl. This also suggests an existence of other 
components than the optical component. 

3.2. CCD chip-misalignment 

Next, we examine PSF eUipticities remaining after sub- 
tracting the optical aberration component. Let us start 
with visual impression. In Figs. 6 and 7, PSF eUipticities 
of observed data (left), the optical aberration model (cen- 
ter), and residual of those two (cobs — Copt-modoi) sue plot- 
ted. Note that the mosaicked data consists of 2-rowx5- 
column CCD chips with masked outer region (beyond 15 
arcmin radius from the field center). In those residual 
plots, there are two apparent features which should be no- 



Fig. 7. Same as Fig. 6 but for another data shown in Fig. 2. 



ticed: one is discontinuities between chips, and the other 
is ripple like patterns. The latter is a characteristic fea- 
ture of PSF ellipticites caused by atmospheric turbulence 
and is also observed in CFHT data (Heymans et al. 2012) 
and in numerical simulation (Jee & Tyson 20011; Chang 
et al. 2012; 2013; also see Appendix 2), which we will 
examine in detail in the next subsection 3.3. 

In this subsection, we focus on PSF eUipticities caused 
by misalignment of CCD chips. To do so, we first need 
to create data that are minimally affected by the optical 
aberration and atmospheric turbulence. The former could 
be minimized by using data with small values of both \as \ 
and a A. The latter could be suppressed by stacking many 
images as the amplitude of the atmospheric PSF decreases 
with the exposure time as (e^)^/^ ex T~^^^ (de Vries et al 
2007; Chang et al. 2012; Heymans et al. 2012; see also 
Appendix 2). Thus we generate a coadded image of 52 
shots of 60-scc exposures taken from G5 group where the 
effect of the optical aberration is in a minimum level and 
a largest number of shots (52 shots) is contained. The 
stacking is done using SWARP, and the object detection 
and ellipticity measurement are done with the same man- 
ner as one described in §2. 

The PSF eUipticities of the coadded data are shown in 
Fig. 8. One can see PSF discontinuities between chips 
even in the observed PSF (top-left panel) as well as in the 
residual PSFs (top-right panel). We attempt to extract 
the chip-misalignment component from residual PSF el- 
lipticities by fitting the residual PSFs with the 2nd order 
bi-polynomial function on the chip-by-chip basis. Derived 
PSF models are shown in the bottom-right panel of Fig. 8 
in which it is seen that characteristic features of chip- 
misalignment PSF eUipticities (such like discontinuities 
between chips and chip-scale smooth variation) are well 
reproduced by the model. We measure the aperture mass 
variances from the PSF eUipticities and show in Fig. 9. 
Note that the B-mode variances are too small to analyze 
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Fig. 8. EUipticity magnitude for each component is shown: the left and right set of panels are for ei and 62, respectively. Results 
from a coaddod data of 52 60-sec exposures from G5 group are shown. For each set of panels, top-left panel shows an observed 
ellipticity, bottom-left panel is for the aberration model, top-right panel shows the residual of those two (eo^s ~ Eopt— model), a-nd 
bottom-right panel shows the chip-misalignment model described in the section 3.2. Note that the color scale of each column is 
shown in top of each column. 

to the optical aberration and thus its amplitude may vary 
with deviation in the focal plane position. Therefore the 
above value should be considered as a characteristic size. 
It is also found from Fig. 9 that the aperture mass vari- 
ance of the residual PSFs (eobs ~ eopt-modoi) is slightly 
larger than that of the chip-misalignment PSF model. The 
amplitude of the excess is roughly {M^p) ^ 10~^. An origin 
of this excess is not clear, it could be due to imperfection 
in the models and/or due to atmospheric PSFs and/or due 
to other causes such like pixelization effects caused in a 
process of resampling (Hamana & Miyazaki 2008). If the 
latter is the case, it can be settled by recently developed 
shape measurements schemes that do not involve a resam- 
pling process (e.g., Miller et al. 2007; 2013; and Miyatake 
et al 2013). Thus we leave it for future work. 

3.3. Atmospheric turbulence 

Let us turn to the PSF ellipticity originated from at- 
mospheric turbulence. Properties of the atmospheric PSF 
ellipticity has been examined experimentally (Wittman 
2005; Heymans et al. 2012) or using numerical simula- 
tions (de Vries et al 2007; Jee & Tyson 20011; Chang 
et al. 2012; 2013). A major characteristic of the atmo- 
spheric PSF ellipticity found from previous studies is that 
the amplitude of PSF ellipticities decreases with the ex- 
posure time as [e?)^!"^ a^T~^l'^ (de Vries et al 2007; Chang 
et al. 2012; Heymans et al. 2012). Also the shape of the 
power spectrum and aperture mass variance of the atmo- 
spheric PSF ellipticities is investigated numerically (see 
Appendix 2). 

We investigate the atmospheric PSFs using data from 
the G5 group. In top-panels of Fig. 10, we show the aper- 
ture mass variances for different exposure times; from top 
to bottom, single 30-sec, 60-sec exposures, coadded 4x60, 
13x60 and 52 x 60-sec exposures. As is shown in the plots. 




0[arcmin] 

Fig. 9. The E-mode aperture mass variance of PSF elliptic- 
ities measured from the coadded data shown in Fig. 8. Solid 
line is for the observed PSFs, dashed line is for the aberra- 
tion model, the dotted line is for the residual of those two 
(eobs — Gopt- model) and long-dashed for chip-misalignment 
PSF model. 

properly and thus we do not present them. Although our 
model may not extract only the chip-misalignment PSFs 
but may be affected by other sources of PSF ellipticity, 
it may be still reasonable to consider that the model al- 
lows us to assess an approximate amplitude of the chip- 
misalignment PSF ellipticies. We read from the Fig. 9 that 
the amplitude of the chip-misalignment PSF ellipticies is 
roughly (M^ ) ^ 10^^ x G. It should be, however, noted 
that the chip-misalignment PSF is a phenomenon related 
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Fig. 10. Top-panels: the aperture mass variances of PSF el- 
lipticities measured from data in G5 group for different ex- 
posure times; from top to bottom, single 30-sec, 60-sec expo- 
sures, coadded 4, 13 and 52 60-sec exposures. Bottom-panels: 
differences between aperture mass variances and one mea- 
sured from the deepest coadded data. 

the amplitude of the aperture mass variance decreases as 
the exposure time increases. However the minimum am- 
phtude of E-mode is set by the optical aberration PSFs 
and chip-misalignment PSFs as shown in the last two sub- 
sections, whereas the B-modc variance decreases down to 
(Mj) ^ 10"^. This indicates that the optical and chip- 
misalignment PSF largely contribute to E-mode PSF el- 
lipticity. 

In order to evaluate the amplitude of the atmospheric 
PSF ellipticities separately from other components, we 
make a working hypothesis that the amplitude of the 
deepest data (52 x 60-sec coadded data) represents an am- 
plitude of (quasi-)static PSF components (such like the 
optical and chip- misalignment PSFs), and the difference 
from that gives an estimate of atmospheric PSFs ellip- 
ticities. In the bottom-panels of Fig. 10, we plot the 
differences where one finds that the aperture mass vari- 
ances have both the exposure time dependence and sharp 
similar to those obtained from the numerical simulation 
(Appendix 2). Also, E- and B-niode variances are in the 
similar amplitude as expected from the numerical simula- 
tion. Although the above results are derived relaying on 
the working hypothesis, those nice agreements would sup- 
port reasonableness of the hypothesis. We may therefore 
say that we successfully capture characteristic features of 
the atmospheric PSF ellipticities in Suprimc-Cam data. 



4. Testing PSF correction scheme using numeri- 
cal simulations 

Here we investigate how well the PSF ellipticities can 
be corrected by the standard correction scheme which we 
describe below. The aim of this study is twofold: One is 
to understand what sets a lower limit of the anisotropic 
PSF correction, and the other is to evaluate an amplitude 
of the residual PSF ellipticities. To do this, we use both 
the dense star field data and mock data that we describe 
below. 

Before going into details, we briefly describe a basic pro- 
cedure of the anisotropic PSF correction implemented in 
actual weak lensing studies so far (see for details, Heymans 
et al 2006; Massey et al 2007 and references therein): 
Firstly, PSF ellipticities (or other estimators of PSF) are 
measured from high-SN star images. In usual weak lensing 
analyses, a typical number density of stars used for PSF 
sampling is ~ 1 arcmin"^. Secondly, the spatial variation 
of PSF ellipticities is modeled by fitting PSF ellipticities 
with an analytic function such like the polynomial. Then 
finally, an artificial shape deformation in galaxy images 
caused by the PSF is corrected using the PSF model deliv- 
ered by the analytic function. Although an actual imple- 
mentation of the above procedure depends on technique of 
weak lensing analysis (Heymans et al 2006; Massey et al 
2007), there is one common limitation which comes from 
a sparse sampling of stars, namely, the spatial variation 
of PSFs on scale smaller than the mean star separation is 
poorly sampled, and as a result, small scale components 
in the PSF anisotropy are hardly corrected. 

Our PSF correction scheme is based on the KSB algo- 
rithm (Kaiser et al 1995; Hamana et al. 2003; Heymans 
et al 2006) and we refer the reader to the above refer- 
ences for details. Here we only describe some details in 
implementation which are specific to this study. We use 
bi-variable polynomial function for modeling the spatial 
variation of PSF ellipticites. We divide data into 2x2 
regions or into each CCD chip, and a PSF model is gen- 
erated for each sub-region. The order of polynomial is set 
4th for 2 X 2-division and 2nd or 3rd for the chip-basis case. 
The 2 X 2-division is taken because in our actual cosmic 
shear analysis the PSF correction is made on a sub-region 
basis with a similar area, whereas the chip-basis analy- 
sis is taken to test its ability of improving the anisotropic 
PSF correction. 

Let us start with a simpler situation using mock PSF cl- 
lipticity catalogs. The mock data have the same geometry 
and number of objects as those of the dense star field data, 
that is 700 "star-role" objects and 30,000 "galaxy-role" 
objects (see §2). Instead of generating mock image, we 
simply generate object catalogs. Each object has items; 
the positions, and 3-component PSF ellipticities including 
optical, atmospheric, and CCD-misalignment components 
denoted by Copt, Catm and Cchip, respectively. The optical 
component is made by using the optical aberration moti- 
vated model, cq. (7), introduced in Appendix 1, and we 
adopted model parameters for G5 group in §3.1. For the 
CCD-misalignment component, we adopt the 2nd-order 
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Fig. 11. E-mode aperture mass variances of PSF ellipticities 
before and after anisotropic PSF correction obtained from 
the mock simulations. Dashed lines are for three input PSF 
models including the optical (black), atmospheric (green) and 
the CCD-misalignment components. The long-dashed lines is 
for total of these three components. The solid line shows 
one after the correction based on the 2 X 2-division, whereas 
the dotted line is for the chip-basis. From left to right, the 
exposure time is 1, 5 and 25 minutes, respectively. 

polynomial models obtained in §3.2. The atmospheric 
component is modeled by the random field with the power- 
law power spectrum of the power-law index of —3 based 
on our finding (see §2 and Appendix 2). Its amplitude is 
set by the empirical relation of eq. (A4) for the exposure 
times of 1, 5, and 25 minutes. We compute the total el- 
lipticity by simply adding the three components, instead 
of treating in the convolution manner. This is a reason- 
able approximation because each cllipticity component is 
small. Then, the anisotropic PSF correction is made for 
mock catalogs in the procedure described above. 

The results from the mock simulation are presented in 
Fig. 11, from which the following three points are seen. 
First, as mentioned above, the mean separation of stars 
(~ 1 arcmin), which are used to model the spatial varia- 
tion of PSFs, sets a fundamental limitation on the PSF 
correction. In the case of the aperture mass, its win- 
dow function is most sensitive to fluctuations on scales 
~ 1/5 of the aperture size (Schneider et al. 1998), This 
defines a characteristic aperture scale of ~ 5 arcmin, and 
on scales larger than that the correction works very well 
but on smaller scale it does very poorly. This is espe- 
cially the case for the atmospheric PSF. Second, in spite 
of the above limitation, the optical components is well 
corrected even on smaller scales. The reason of this is 
the smoothness of the optical component; as has been 
found in §3.1, the spatial variation of the optical PSFs 
can be well fitted by a polynomial-based model. Thus 



Fig. 12. E- (left) and B-mode (right) aperture mass vari- 
ances of PSF ellipticities before (dashed curves) and after 
(solid curves) the anisotropic PSF correction. The results ob- 
tained from the G5 group arc presented for various exposure 
time. The PSF correction is made on a basis of 2 X 2 sub-re- 
gions except for the magenta curve which is for the chip-basis 
correction. 

even with a sparse sampling, it can be well modeled with 
the polynomial function. Third, for 2 x 2-division case, the 
CCD-misalignment component combined with the optical 
component sets a limit on the anisotropic PSF correction. 
This is a natural consequence that even each individual 
component is smooth and is modeled by polynomial func- 
tion, the mixed PSFs may be no longer simple and thus 
may not be corrected by the polynomial model. This can 
be partly avoided by adopting the chip-basis correction as 
is seen in the right panel of the Fig 11, though the resid- 
ual variance does not reach to the level of the atmospheric 
PSF. A further improvement of the PSF correction may 
be achieved by recently proposed correction schemes (e.g., 
Chang et al. 2012; Bergc et al. 2012; Gentile, Courbin & 
Meylan 2013), and we leave it for future study. The above 
three findings provide us with a clue to interpret a result 
obtained from an analysis of real data. 

Let us turn to the analysis of the dense star field data. 
In Fig. 12, the aperture mass variances measured from the 
G5 group data (defined in Fig. 4) are presented for differ- 
ent exposure times, where the anisotropic PSF correction 
is made on a basis of 2 x 2 sub-regions except for the ma- 
genta curve to which the chip-basis correction is applied. 
From the figure, the following four findings are obtained: 
Firstly, it is found that after the PSF correction E- and 
B-mode are almost equally partitioned, regardless of prop- 
erties of raw PSF ellipticities. A plausible reason of this 
is that E- and B-mode are mixed up in the process of the 
PSF correction. Secondly, the result of 60-sec exposure 
is in a good agreement with that of the mock simulation. 
But at a closer look at larger scales, it is found that the 
residual aperture mass variance falls less quickly than the 
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mock result, indicating tliat the spatial variation of real 
PSF ellipticities is more complex than one assumed in the 
mock simulation. Thirdly, as expected from the result of 
the mock simulation, the residual aperture mass variance 
decreases with increasing the exposure time. However, the 
decrement is less than what expected from the mock sim- 
ulation. Finally, it is found from the comparison of two 
different correction schemes for 52-coadded data that the 
chip-basis correction little makes an improvement over the 
other case. It is speculated from the last two points com- 
bined with the findings from the mock simulation (that 
is, an existence of a non-smooth component sets a lower 
limit of the correction on scales smaller than mean separa- 
tion of stars) that there exists unidentified PSF ellipticity 
component (s) whose amplitude is greater than the chip 
component found in §3.2 (see also Fig. 9). Although an 
origin of the unknown PSF component is not clear, the 
above result provides us with an empirical estimate of the 
"best performance" of the PSF correction scheme that we 
adopted in this paper. 

5. Residual PSF anisotropies in Suprime-Cam 
cosmic shear data 

Here wc examine effects of the residual PSF anisotropies 
on Suprime-Cam cosmic shear data. To do so, we measure 
the E- and B-mode shear aperture mass variance from 5.6- 
degree^ Suprime-Cam deep imaging data, and compare 
the B-mode shear variance with the residual PSF elliptic- 
ities measured from the dense star field in §3. 

5.1. Suprime-Cam data and weak lensing analysis 

Wc use the same data as those used in Hamana et 
al (2012) in which basics of data and data analyses are 
described. Therefore here we only describe points spe- 
cific to this work. Suprime-Cam i'-band data are cor- 
rected from the data archive, SMOKA, under the follow- 
ing three conditions: data are contiguous with at least 
four pointings, the total exposure time for each point- 
ing is longer than 1800 sec, and the seeing full width at 
half-maximum (FWHM) is better than 0.65 arcsec. Four 
data sets (named by SXDS, COSMOS, Lockman-hole and 
ELAIS-Nl) meet these requirements. The exposure times 
of individual exposures range from 120 and 420 seconds. 
Data reduction, mosaic stacking and object detection were 
done with the same procedure described in §2. The effec- 
tive area after masking regions affected by bright stars is 
5.57 degree^. Depth of the coaddcd images varies among 
four fields, and it is found that the number counts of faint 
galaxies are saturated at AB magnitude of i' = 25.2 — 25.5. 

For weak lensing measurements, we adopt the so-called 
KSB method described in Kaiser et al. (2005), Luppino 
& Kaiser, and Hoekstra et al (1998) with some modifica- 
tions being made following recent developments (Heymans 
et al. 2006a; Massey et al. 2007). Stars are selected in a 
standard way by identifying the appropriate branch in the 
magnitude half-light radius (r/^) plane, along with the de- 
tection significance cut S/N > 10. The number density of 
stars is found to be ~ 1 arcmin"^ for the four fields. We 
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Fig. 13. Same as Fig. 4 but for data used for tfic cosmic sliear 
analysis. The optical model fitting is applied to individual 
exposures. The exposure time of each exposure is 120 — 420 
seconds. Different symbols are for different fields. 

only use galaxies met the following three conditions, (i) 
the detection significance of S/N > 3 and nu > 10 where 
nu is an estimate of the peak significance given by hfind- 
peaks, (ii) r^ is larger than the stellar branch, and (iii) the 
AB magnitude is in the range of 22 < i' < 25. The number 
density of resulting galaxy catalog is '^ 23.5 arcmin"^ and 
is quite uniform among four fields. We measure the shape 
parameters of objects by getshapes of IMCAT. The PSF 
correction is done on a sub-field basis, where the coadded 
data were divided into sub-fields of about 15 x 15 arcmin^ 
(approximately one fourth of the Suprime-Cam's field-of- 
view). The spatial variation of PSF ellipticites measured 
from stars are modeled with the 4th order bi-polynomial 
function, and the PSF model combined with the smear 
polarizability tensor (Psm) of each galaxy is used to cor- 
rect the PSF ellipticities of galaxies. In KSB formalism, 
the shear (we denote by 7) is related to the observed ellip- 
ticity through the shear polarizability tensor, P^, which 
is evaluated by the smoothing and weighting method de- 
veloped by Van Waerbeke et al. (2000; see Hamana et al 
2003 for an actual implementation of this method) . 

Before presenting the shear aperture mass variance 
measured from coadded data, it is worth looking into PSF 
ellipticities of individual exposures. To do so, we gen- 
erate mosaicked data of individual exposures and select 
stars in the same procedures as described above. The 
spatial variation of PSF ellipticities measured from the 
stars are fitted to the optical model, eq. (7), and the am- 
plitude of the symmetric and asymmetric components of 
the model are evaluated (see §3.1 and Appendix 1 for de- 
tails). The results are plotted in Fig. 13, in which one 
may see "V-shaped" trend similar to the results obtained 
from the dense star field data (Fig. 4), though the scat- 
ter is somewhat larger for this case. This similarity im- 
plies that knowledge of the PSF anisotropy learned from 
the dense star field analysis (§3) is applicable to cosmic 
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Fig. 14. Filled circles show E- (left) and B-mode (right) 
shear aperture mass variances measured from faint galax- 
ies. The error bars present the statistical error computed 
from 100 randomized realizations. The crosses are for un- 
corrected- shears (the galaxy shears to which the correction 
for the PSF anisotropy is not applied). The dashed curves 
shows variances of the residual PSF ellipticities (but scaled to 
the shear by multiplying a factor ({Psm / P-y) / (Pgrn))^ — 1-5^) 
measured from the G5 group of the dense star field data (sec 
§3 and Fig. 12); the upper and lower curves are for 13 X 63 and 
52 X 60-scc exposure, respectively. The dotted curve shows the 
ACDM prediction based on the 7-year WMAP cosmological 
model with the source galaxy redshift distribution inferred 
from COSMOS photometric redshift catalog. 

shear analyses. It should be however noticed that there 
is one component that is not taken into consideration in 
the dense star field analysis, namely stacking of multiple 
dithered exposures. It is obvious that stacking of dithered 
exposures makes the spatial variation of PSFs complicated 
which may result in a poorer anisotropic PSF correction. 
We will return to this point later. 

5.2. Cosmic .shear aperture mass variance 

We compute the E- and B-mode aperture mass vari- 
ance of galaxy shears through the two-point correla- 
tion functions using cqs. (3)-(5) with the summations 
in eq. (5) being replaced with the weighted summa- 
tion, ^WiWj^i'^j/^WiWj where Wi is the weight for i- 
th galaxy. In addition to the shear, we compute the 
uncorrected- shear, that is, the galaxy shear to which the 
correction for the PSF anisotropy is not applied. The re- 
sults are presented in Fig. 14. The error bars indicate the 
root-mean-square among 100 randomized realizations, in 
which the orientations of galaxies are randomized, and 
presumably represent the statistical error coming from 
the galaxy shape noise. It is found from Fig. 14 that 
the B-mode shear variance is non-zero on the aperture 
scale of ^ 10 arcmin. It should be noted that adja- 
cent bins are correlated. In order to compare the shear 
variances with those of the residual PSF ellipticites, we 
transform the star ellipticity to the galaxy shear by mul- 
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Fig. 15. Redshift distribution of galaxies used for the cosmic 
shear analysis derived adopting the COSMOS photometric 
redshift catalog (Ilbert ct al 2009). The solid curve shows 
the fitting model by eq. 8 with the parameters denoted in the 
panel. 

tiplying a factor {Psm / P-y) / (Psm) which is found to be 
~ 1.5 (thus 7* ~ 1.5e*). The results from the G5 group 
of the dense star field data with the 2x2 sub-region basis 
correction (see §3 and Fig. 12) are also shown in Fig. 14. 
The upper and lower curves arc for 13 x 60 seconds and 
52 X 60 seconds exposure, respectively. Since the total ex- 
posure times of four cosmic shear fields are 2400 — 6000 
seconds, the lower curve could be considered as a "6est per- 
formance^^ of our PSF ellipticity correction scheme. The 
B-mode shear signals at smaller scales {6 ^ 10 arcmin) are 
larger than this expectation. A possible reason of this ex- 
cess is that stacking of dithered multiple exposures which 
is not involved in the dense star field analysis. This issue 
will be investigated in detail in a future study in which 
we will adopt a PSF correction scheme on an individual 
exposure basis. On the other hand, the B-mode variance 
drops at larger scales (0 ^ 15 arcmin), which is in nice 
agreement with one expected from the dense star field 
analysis. It is however noticed that in the case of the dense 
star field data, the turnaround aperture scale is about 5 
arcmin which is naturally set by the mean separation of 
stars (~ 1 arcmin) (see §4), whereas it is about 10 arcmin 
for the cosmic shear data. The reason of this is unclear, 
it could be related to an excess B-modc uncorrected shear 
variance observed at '--^ 10 arcmin. 

Since the anisotropic PSF correction works well on 
larger scales, it is interesting to compare the E-mode vari- 
ance with the theoretical prediction on those scales, we 
derive the redshift distribution of source galaxies, which is 
an essential ingredient for computing a theoretical model, 
in the following manner: We utilize the COSMOS field 
data. We merge the galaxy catalog used for the cosmic 
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shear measurement witli the COSMOS photometric red- 
shift catalog (Ilbert et al 2009), and compute the redshift 
distribution adopting the photometric redshift, which is 
presented in Fig. 15. We fit the redshift distribution with 
the following function (Fu et al. 2008) 



n{z) ^ A 






(8) 



where the normalization A is determined by imposing 
J dz n{z) ~ 1 within the integration range < z < 6. We 
find that a parameter set {a,h,c) = (0.7,6.0.7) gives a rea- 
sonably good fit to the data as shown in Fig. 15. The 
mean redshift of this model is 1.07 which is in a very 
good agreement with the value derived from the photo- 
metric redshift data {{zp) = 1.06). We plot in Fig. 14 the 
ACDM prediction based on the 7-year WMAP cosmologi- 
cal model (Komatsu et al. 2011). For the nonlinear CDM 
power spectrum, we use the HALOFIT model by Smith et 
al. (2003). It is seen in Fig. 14 that the measured E-mode 
shear variance is in a nice agreement with the theoretical 
model on the scales where the anisotropic PSF correction 
works well (0 > 15 arcmin). One may be interested in a 
constraint on the cosmological parameters to be placed 
from the data. However, to derive that, there are some is- 
sues to be carefully taken into consideration; for example, 
an effect of the noise bias (Refregier et al. 2012; Melchior 
& Viola 2002; Okura & Futamase 2012) and the covari- 
ance matrix of cosmic shear correlation function includ- 
ing non-Gaussian contributions (Sato et al., 2009; 2011). 
Those are beyond of the scope of this paper, and will be 
examined in a future work. 

6. Summary and discussion 

We have analyzed the anisotropic PSF of Suprime- 
Cam data utilizing the dense star field data. We have 
decomposed the PSF ellipticities into three components, 
the optical aberration, atmospheric turbulence and chip- 
misalignment in an empirical manner, and have assessed 
the amplitude of each component. Then, we have tested 
a standard method for correcting the PSF ellipticities 
against mock simulation which is based on models of the 
PSF ellipticities obtained from the dense star field data 
analysis. Finally we have examined the impact of the 
residual PSF ellipticities on cosmic shear measurement. 

Wc have found that the optical aberration has the 
largest contribution to the PSF ellipticities of long- 
exposure (say, longer than 10 minutes) data. We have also 
found that the spatial variation of the optical PSFs can 
be modeled by the simple three-component model which 
is based on the lowest-order aberration theory described 
in Appendix 1. It is also found that the optical PSF ellip- 
ticities vary smoothly on the focal plane, and thus are well 
corrected by the standard correction method in which the 
spatial variation is modeled by a polynomial function. 

The PSF ellipticities after subtracting the optical PSF 
model ellipticities show discontinuities between CCD 
chips, which indicates that the misalignment of chips on 
the focal plane induces an additional PSF ellipticity. In 



order to make a crude estimation of the amplitude of chip- 
misalignment PSF ellipticities, we have fitted the residual 
PSF ellipticities to the 2nd order bi-polynomial function 
on chip-by-chip basis. It turns out that the amplitude of 
the chip-misalignment PSF ellipticities is lower than the 
optical component, though it is not negligible. It is found 
from the mock simulation that the chip-misalignment 
component combined with the optical component puts a 
lower limit on the capability of the PSF correction, which 
can be, in principle, avoided by employing a chip-basis 
correction scheme. 

We have investigated properties of PSF ellipticities orig- 
inated from the atmospheric turbulence using the numer- 
ical simulation of wave propagation through atmospheric 
turbulence under a typical weather condition of Subaru 
telescope at Mauna-kea (Appendix 2). From the simula- 
tion results, we have evaluated the power spectrum and 
aperture mass variance of atmospheric PSF ellipticities, 
and have derived power-law fitting functions of them. As 
was already pointed out in literatures (Wittman 2005; de 
Vries et al 2007; Jee & Tyson 20011; Chang et al. 2012; 
2013 Heymans et al. 2012), the RMS amplitude of the at- 

— 1/2 

mospheric PSF ellipticities decreases as Texp ■ We com- 
puted the aperture mass variance of the dense star field 
data for various exposure time (from 30-sec to 52 x 60- 
sec). The amplitude of the atmospheric PSF ellipticities 
was evaluated by assuming the deepest data represents an 
amplitude of (quasi-)static PSF components, and the dif- 
ference from that is due to the atmospheric PSFs elliptici- 
ties. The results are found to be in a reasonable agreement 
with the simulation results. Therefore, we may conclude 
that the aperture mass variance of the atmospheric PSF 
ellipticies for a long-exposure data (say, longer than 10 
minutes) , is at least one order of magnitude smaller than 
that of the optical PSF ellipticities. Since the atmospheric 
PSFs are not smoothly varying component, its correction 
suffers from the common limitation that the PSF correc- 
tion on scales smaller than the mean star separation works 
very poorly because of the spatial variation of PSFs on 
those scales is hardly sampled. 

The above findings provide us a clue to develop an op- 
timal PSF interpolation scheme: It is found that the spa- 
tial variation of PSF ellipticities consists of two compo- 
nents; one is a smooth and parameterizable component 
arisen from the optical aberration and chip- misalignment, 
and the other is a non-smooth and stochastic component 
arisen from the atmospheric PSFs. The former can be 
modeled with a parametric model as shown in this pa- 
per. Also it has been argued that an interpolation scheme 
based on the principle component analysis is effective for 
such a case (Jarvis & Jain 2004; Jee & Tyson 2011; see 
also Miyatake et al 2013 and Lupton et al 2001 for an 
actual implementation for Suprime-Cam data reduction 
pipeline). On the other hand, it has been shown in Berge, 
et al (2012) and Gentile et al. (2013) that a local type 
correction scheme such as the radial basis functions and 
Kriging works well for the atmospheric PSFs. Apparently, 
a hybrid interpolation scheme, in which the above two 
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types ol interpolation schemes are optimally incorporated, 
is a strong candidate for achieving a better PSF correc- 
tion. 

We have examined effects of the residual PSF 
anisotropics on Suprime-Cam cosmic shear data. We have 
compared the B-mode shear variance measured from 5.6- 
degree^ I'-band data with the residual PSF ellipticities 
(but being properly transformed into shear) measured 
from the dense star field data, which can be considered as 
the "best performance" of our PSF ellipticity correction 
scheme. It is found that the shape and amplitude of B- 
mode shear variance are broadly consistent with those of 
the residual PSF ellipticities. This indicates that most of 
the sources of residual systematic are understood, which is 
an important step for cosmic shear statistics to be a prac- 
tical tool of the precision cosmology. However it is also 
found that the B-mode shear amplitude at scales ^ 10 
arcmin are systematically larger than the residual PSF. 
A reason of this excess is unclear, one possible reason is 
that stacking of dithered multiple exposures which is not 
involved in the dense star field analysis. Such stacking 
related issues may be avoided by employing a weak lens- 
ing shape measurement scheme on the basis of individual 
exposures (e.g.. Miller et al. 2007; 2013; and Miyatake 
et al 2013), which combined with a hybrid interpolation 
scheme will be addressed in a future work. 

We greatly thank Y. Miller for many fruitful discus- 
sions on various aspects of this work. We thank M. 
Takada, M. Oguri and H. Miyatake for useful discus- 
sions. We thank Y. Utsumi for help with photometric 
calibration of Suprime-Cam data. We also thank L. Van 
Waerbekc and R. Mandelbaum for valuable comments on 
on an earlier manuscript which improved the paper. We 
would like to thank M. Britton for making the Arroyo 
software available. Numerical computations in this pa- 
per were in part carried out on the general-purpose PC 
farm at Center for Computational Astrophysics, CfCA, of 
National Astronomical Observatory of Japan. This work 
is based in part on data collected at Subaru Telescope 
and obtained from the SMOKA, which is operated 
by the Astronomy Data Center, National Astronomical 
Observatory of Japan. This work is supported in part 
by Grant-in-Aid for Scientific Research from the JSPS 
Promotion of Science (23540324). 

References 

Bartelmann M., Schneider P., 2001, Physics Report, 340, 291 
Berge J., Price S., Amara A., Rhodes J., 2012, MNRAS, 419, 

2356 
Bertin E., Arnouts S., 1996: A&AS, 317, 393 
Bertin E., Mellier Y., Radovich M., Missonnier G., Didelon P., 

Morin B., 2002, ASP Conference Proceedings, 281, 228 
Bertin E., 2006, ASP Conference Series, 351, 122 
Bridle S. et al. 2010, MNRAS, 405, 2044 
Britton, M. C, 2004, Proc. SPIE, 5497, 290 
Catelan P., Kamionkowski M., Blandford R. 2001, MNRAS, 

320, 7 
Chang et al, 2012, MNRAS, 427, 2572 



Chang et al, 2013, MNRAS, 428, 2695 

Crittenden R., Natarajan P., Pen U., Theuns T. 2002, ApJ, 

559, 552 
Croft R., Metzler C. 2000, ApJ, 545, 561 
de Vries W. H., Olivier S. S., Asztalos S. J., Rosenberg L. J., 

Baker K. L., 2007, ApJ, 662, 744 
EUerbroek B. L., Rigaut F. J., 2000, Proc. SPIE, 4007, 1088 
Fu et al., 2008, A&A, 479, 9 

Gentile M. Courbin P., Meylan G., 2013, A&A, 549, Al 
Hamana T., et al, 2003, ApJ, 597, 98 
Hamana T., Miyazaki S., 2008, PASJ, 60, 1363 
Hamana T., Oguri M., Shirasaki M., Sato M., 2012 MNRAS 

425, 2287 
Hardy J. W., 1998, "Adoptive Optics for Astronomical 

Telescopes" Oxford University Press 
Heymans C, et al. 2006, MNRAS, 368, 1323 
Heymans C, White M., Heavens A., Vale C, van Waerbeke 

L., 2006, MNRAS, 371, 750 
Heymans C, Rowe B., Hoekstra H., Miller L., Erben T., 

Kitching T., van Waerbeke L., MNRAS, 421, 381 
Hilbert S., Hartlap J., White S. D. M., Schneider P., A&A, 

499, 31 
Hirata C, Seljak U., 2004, Phys. Rev. D, 70, 063526 
Hoekstra H., 2004, MNRAS, 347, 1337 
Hoekstra H., Franx M., Kuijken K., Squires G., 1998, ApJ, 

504, 636 
Hoekstra H., Jain B., 2008, Annual Review of Nuclear and 

Particle Systems, 58, 99 
Huterer D., Takada M., Bernstein G., Jain, B., 2006, MNRAS, 

366, 101 
Ilbert O., et al., 2009, ApJ, 690, 1236 
Jarvis M., Jain B., 2004, submitted to ApJ (arXiv:astro- 

ph/0412234) 
Jarvis M., Schechter P., Jain B., 2008, submitted to PASP 

(arXiv:0810.0027) 
Jee M. J., Tyson J. A., 2011, PASP, 123, 596 
Jee, M J., Blakeslee J. P., Sirianni M., Martel A. R., White R. 

L., Ford H. C, 2007, PASP, 119, 1403 
Jing Y., 2002, MNRAS, 335, 89 

Kaiser N., Squires G., Broadhurst T., 1995, ApJ, 449, 460 
Kaiser N., Torny J. L., Luppino, G., A., 2000, PASP, 122, 768 
Kitching T. D., et al, 2012a, MNRAS, 423, 3163 
Kitching T. D., et al, 2012b, ApJS submitted 

(arXiv/2012.1979) 
Komatsu E., et al., 2011, ApJS, 192, 18 
Luppino G. A., Kaiser N., 1997, ApJ, 475, 20 
Lupton R., Gunn J. E., Ivezic Z., Knapp G. R., Kent S., 

Yasuda N., 2001, in Harnden F. R. Jr, Primini F. A., Payne 

H. E., eds, ASP Conf. Ser. Vol. 238, Astronomical Data 

Analysis Software and Systems X. Astron. Soc. Pac, San 

Francisco, p. 269 
Massey, R., et al., 2007, MNRAS, 376, 13 
Mellier Y., 1999, Annu. Rev. Astron. Astrophys., 37, 127 
Melchior P., Viola M., 2012, MNRAS, 424, 2757 
Miller L., Kitching T. D., Heymans C, Heavens A. F., van 

Waerbeke L., 2007, MNRAS, 382, 315 
Miller L., et al, 2013, MNRAS, 429, 2858 
Miyatake, H., et al, 2013, MNRAS, 429, 3627 
Miyazaki S. et al., 2002, PASJ, 54, 833 
Munshi D., Valageas P., van Waerbeke L., Heavens A., 2008, 

Physics Report, 462, 67 
Okura Y., Futamase T., 2012, subitted to ApJ 

(arXiv:1208.3564) 
Ouchi, M., et al. 2004, ApJ, 611, 660 



No. 



Anisotropic PSF of Suprime-Cam 



13 



r 


- 


f> 


& 


. 


A 


<! 




-? 


: ^ 


^ 


E 


fr 


»; 


A 


^ 


^ 


<? 


\', 


Vj 


© 


6 


6 


& 


(? 


■;^ 


rS 


■^J 


-o 


^ 


, 


• 


. 


137 


e? 


(^ - 


- «© 


-c 


o 


- 




• 


O 


(> 


c> - 


A. 


.fi 


^ 


e 


. 


, 


.» 


Ci 


("* - 


» 


i 


^ 


f 


« 


« 


•4 


X 


i 


^ 


ii 


p 


v> 


V 


•« 


■ij 


~4 


i 


i 


i 


;» 


V 


• 


V 


1^ 


'■^ 


4 _ 



4 « 






,u ,i' s r I 1 4 !si 



4 


-t» 


^ 


!^ 


9 


^' 


j? 


P' 


£>: 


_ >.«, 


h 


Q 


<? 


9 


P 


P 


fi) 


yC,.' _ 


^ 


A 


!^ 


•i 


« 


p 


e 


^) 


Ji-\ 


- (^ 


©^ 


Ci 


« 


• 


» 


^ 


-© 


4) - 


-o 


(> 


<» 


- 




• 


« 


<> 


<)- 


- -V. 


-0- 


(17 


• 


• 


. 


^ 


■© 


«)- 


- >?■ 


O' 


(? 


^ 


s 


6 


o 


^;, 


■■-■■ 


: -^ 


■ff 


-9 


« 


6 


ki 


b 


ft 


v„ ; 


" ^ 


1 


1 


A 


A 


A 


^ 


b 


r " 



I? ^ 

* 

« ' 

« • 

« > 

% » 

% % 



• » % 

• • i 

• t I 



^ 





















J 


; ' 


' 


■1 


^ 


- 


-^ 


l" 


" 


^ 




" 


'• 


'' 


' 


" 


r 


•» 


(? 


' 




" 


-T 


c 


' 


' 


." 


C 


i' 


' 




"? 


t 


' 




■^ 


rr 


' 


' 


f 




- - 


- 


-" 


.' 


7^ 





J^ 





j^ 




,- 


- 


' 


:< 


rf 


(^ 


- 


.' 


-' 




' 


• 


1"' 


' 


i"T 




::? 


' 


r 




- 


^ 


- 


' 


i:? 


' 


' 


(5 


• 




r^ 


' 


' 


" 


1-? 


7 


• 


' 


:5 





/ / / / / 
I / / / 



\ 1 I 

\ \ > 

h \ \ ^ 

\ \ \ 

\ \ \ 



/ / / 







\ \ \ \ 
\ \ \ \ 
\ \ \ \ 



\ \ "^ 



Fig. 16. To-p two panels: spot diagrams for the Seidel coma 
with opposite signs of Wiai. In both cases, the major axis of 
the cllipticity of the PSF are radially oriented with a scaling 
of |e| DC r^. Bottom two panels: spot diagrams for a combi- 
nation of astigmatism and field curvature: Bottom-left panel 
is for a case for \W222 +^^220 1 > |M^22ol where are the bot- 
tom-right panel is for \W222 + W220I < 1^220 1- In both cases, 
the cllipticity parameter is scaled as |e| oc r*. 
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Appendix 1. PSF ellipticity originated from the 
third-order optical aberrations 

We summarize the properties of ellipticity in the PSF 
caused by the third-order optical aberrations with and 
without rotational symmetry. We consider a system with 
a circular pupil. We follow the notation by Thompson 
(2005; see Fig 1 of his paper): H represents the position 
in the image field and without loss of generality we can 
choose {Hx,Hy) — (0,/i). p represents the position in the W ocap cos( 



Fig. 17. Top two panels: spot diagrams for the misalign- 
ment coma with different misalignment directions. Bottom 
two panels: spot diagrams for the misalignment astigmatism 
with different misalignment directions. In both cases, the el- 
lipticity parameter is scaled as |e| ocr^. 

exit pupil with {px,Py) = {psvii(f>, p cos 4>). 

Let us start with the rotationally symmetric system. 
In this case, the wavefront distribution is given by 
(Thompson 2005) 

W = Wo40P'* + Wizihp^ COS<j) + W222h?P^ cos^ <t> 



+W22oh^P^ + W^iih^pcosc, 



(Al) 



Among those five Seidel aberrations, ones that make an 
ellipticity parameter e non-zero are the Seidel coma (the 
second term) and astigmatism (the third term). The 
Seidel coma results in an cllipticity whose major axis 



oriented toward radial direction with lei 



where 



is the distance from the center of image field. On the 
other hand, a combination of the astigmatism and field 
curvature (the fourth term) results in radially elongated 
(if \W222 + W22o\ > IVF220I) or tangentially elongated (if 
IW222 + W220I < |W^22o|) ellipticity parameter with |e| ocr* 
(see Fig. 16). 

Next, we consider the case without rotational symme- 
try, allowing for a tilt and decenter between the axes of 
the exit pupil and image field. To compute the wavefront 
distribution in a perturbative approach, it is customary to 
introduce a vector a with {(Jx,(yy) = (cr sin 0,0" cos 0) which 
represents the decentration of the center of the aberration 
field W with respect to the unperturbed field center. In 
this case, optical aberrations result in a PSF with non- 
zero ellipticity parameter are the misalignment coma and 
astigmatism. 

The wavefront distribution for the misalignment coma 
is given by 

9). (A2) 
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In this case, the PSF has the same coma shape as the 
Seidel coma. However, the PSF does not depends on the 
position in the image field as the above expression does 
not include h, but the PSF (and thus the cUipticity) is 
constant over the image field as shown in the two top 
panels of Fig. 17. 

The wavefront distribution for the misalignment astig- 
matism is given by 

W oc ahp^ cos^{ 



>-e/2). 



(A3) 



In this case, the ellipticity of PSFs depends not only on 
the distance from the center of the image field but also on 
the azimuth angle (we shall denote it by ^), and is written 
by {ei,e2) <x {r'^ cos{iIj — 9),r'^ sm{i/; — 6)) (see Fig. 17). 

Although an actual optical aberration may result in 
more complicated PSF than that expected from the third- 
order optical aberrations, it may be reasonable to build a 
model of the PSF ellipticity caused by optical aberration 
based on the above results. We consider the following form 
of a model (already introduced in eq. (7) but rewritten for 
self-contained) that consists three terms, namely the con- 
stant, axis-symmetric and asymmetric components; 



ei 

62 



Cl 
C2 



E^ 



ris 



cos 2';/' 
sin2V' 



Z^""""^ \ sin(V'-6') 



We determine the model parameters s„^ , a„^ and 9 by the 
standard least square method with observed PSF elliptic- 
ities. In doing so, we take rig = {1,2,3,4} and Ua = {1,2}. 
As was shown in §3, this simple model succeeds very well 
in fitting the global pattern of PSF ellipticities in Suprime- 
Cam data. 

Appendix 2. A numerical study of PSF elliptici- 
ties caused by atmospheric turbulence 

We examine statistical properties of PSF ellipticities 
caused by atmospheric truculcncc using a numerical sim- 
ulation. The aim of this study is to investigate the shape 
and amplitude of the power spectrum of atmospheric PSF 
ellipticities, which is an essential knowledge for under- 
standing statistical properties of PSF ellipticities. Our 
approach is similar to one taken by Jee & Tyson (2011), 
that is to utilize a numerical simulation of wave propaga- 
tion through atmospheric turbulence for evaluating atmo- 
spheric PSFs. 

Here we briefly describe our numerical simulation. Since 
we heavily use a publicity open software Arroyo (Britton 
2004), we refer the reader to the above reference and 
software documents^ for details of models, computational 
methods and their implementation. The atmospheric tur- 
bulence is modeled by eleven-layer frozen Kolmogorov 
screens with the atmospheric model at Mauna-Kea pro- 
posed by EUerbroek & Rigaut (2000). The strength of 
the total atmospheric turbulence is set by the Fried length 




X [arcmin] 



[arcmin] 



X [arcmin] 



http://eraserhead.caltech.edu/arroyo/arroyo.html 



Fig. 18. Ellipticity magnitude obtained from the numerical 
simulation described in Appendix 2: upper and lower panels 
arc for ei and 62, respectively, and the color scales are dis- 
played on the top of those panels. The exposure times are, 
from left to right, 1-sec, 10-sec and 60-sec, respectively. 

To for which we adopt a value proposed by EUerbroek 
& Rigaut (2000); ro = 0.23m at a wavelength of 0.5/im. 
Neither a large scale nor small scale cutoff (the, so-called, 
outer and inner scale) on the Kolmogorov power spec- 
trum is imposed. The wind model given by eq. (3.20) 
of Hardy (1998) is adopted with random wind directions. 
Having the atmospheric model ready, an atmospherically 
disturbed phase function of a wavefront is computed. 
We consider a monochromatic electromagnetic wave with 
wavelength of 0.8//m. Then, an instantaneous PSF is ob- 
tained by Fourier transforming the phase function within 
a pupil for which we assume a circular pupil with 8.2m 
diameter for simplicity. All the above models and compu- 
tations are implemented in Arroyo subroutines. Finally, 
a sequence of instantaneous PSFs arc added to obtain a 
long-exposure PSF. We compute PSFs of 1, 10 and 60 
seconds exposures on 64 x 64 regular grids over a 30' x 30' 
focal plane. 

We generate 20 realizations. The mean FWHM of PSFs 
among the realizations is 0.65 arcsec for 60 seconds ex- 
posure which is in a good agreement with the observed 
value at Subaru telescope (Miyazaki et al. 2002). For 
an illustrative example, one realization is shown in Figure 
18, where the ripple like feature appears. Similar fea- 
tures are observed in CFHT data (Heymans et al. 2012) 
and Suprime-Cam data (see Figs 6 and 7). It should be 
however noted that the appearance varies widely between 
realizations depending on turbulence and wind properties. 

We measure the power spectrum and aperture mass 
variance from the simulation data and present in Figure 
19 and 20, respectively. The following three major find- 
ings are derived from the Figures: First, the amplitude of 
the power spectrum decreases as Pg oc T^^ as was already 
pointed out in previous studies (Wittman 2005; de Vries 
et al 2007; Jee & Tyson 20011; Chang et al. 2012; 2013 
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Fig. 19. The power spectrum of PSF ellipticities measured 
from the simulation results. The dashed, long-dashed and 
solid lines are for 1-sec, 10-sec and 60-sec exposures, respec- 
tively. The error bars, plotted only for the 60-sec result for 
clarify, show the RMS among 20 realizations. The dotted line 
shows the power-law model of Pe(0 = 0-5 ^ l~^ which turns 
out to be a good approximation of the 60-sec result. 




9[arcmin] 



Fig. 20. The aperture mass variance of PSF ellipticities mea- 
sured from the simulation results. The dashed, long-dashed 
and solid lines are for 1-sec, 10-scc and 60-sec exposures, re- 
spectively. The error bars, plotted only for the 60-sec result 
for clarify, show the RMS among 20 realizations. 



Heymans et al. 2012). Second, the power spectrum is ap- 
proximated by a power-law model with the power index 
of —3. A plausible reason of the power-law shape is our 
assumption of the power-law turbulence power spectrum 
without a cutoff. In reality, a large-scale cutoff (the, so- 
called, outer scale) at around 10— 100m is expected (e.g.. 
Hardy 1998), which may result in a cutoff in the PSF ellip- 
ticity power spectrum at a scale around 1 degree. Third, 
the atmospheric turbulence results in E/B-mode almost 
equally partitioned PSF ellipticities. At a closer look at 
Fig. 20, one may however find that the E-mode ampli- 
tude is slightly but systematically larger than B-mode. A 
reason of this is unclear and we leave it for future work. 

In addition to the above findings, we also find 
from Figs. 19 and 20, that models Pe{l) ~ O.S^-^ and 
{Map,±}{d) -- 4 X 10"'^ 61 give a good fit to the 60-sec resuh. 
In what follows, we discuss how those correlation ampli- 
tude depends on other observational parameters includ- 
ing the exposure time and the atmospheric seeing FWHM 
(^atm). Consider a ease D :s> ro where D denotes a diam- 
eter of a telescope, and a monochromatic radiation with 
a wavelength A from a point source. A wavefront from 
a point source is disturbed by atmospheric turbulence, 
and we consider, for simple approximation, the disturbed 
wavefront as segmented rp-sized patches of constant phase 
and random phases between the individual patches (Saha 
2007, and see e.g.. Fig 5.3 of that textbook for an illus- 
tration). The number of patches within a pupil is approx- 
imately Npatch ~ {D/r^f'. A wavefront from each patch 
results in a sharp speckle PSF on a focal plane, and a 
total PSF is viewed as a superposition of randomly dis- 
tributed speckles over an extent of the atmospheric seeing 
size 0atm ^ A/ro (Hardy 1998; Saha 2000; see for illustra- 
tive examples. Fig 3 of Kaiser. Tonry & Luppino 2000; 
and Fig 5 of Jee & Tyson 2011). In this picture, the PSF 
anisotropy is understood as a natural consequence of a 
finite number of speckles. Since the disturbed wavefront 
changes quickly, more different speckles are accumulated 
as Nspec fx Tea;p, and thus the PSF becomes rounder. It is 
found from the above numerical simulation that the RMS 
of PSF ellipticities decreases as {e^y^^ oc Texp , which 
combined with the above mentioned relationships leads 
to the following relationships; 



cxiv-i/2 r-l/2^-1/2 

^^^^spec ^^^exp -^^ patch 



(e2)i/2 

« T^xl^'^1'0 oc Te^yV^atm. 

We test this relationship against numerical simulations 
with different values of ro, from which the above depen- 
dences of tq and 6'atm on (e^)^/^ are verified. Combined 
the last relationship with the power-law fitting models, we 
have 

-1 / /I \ -2 



Pe{l)'^ 0.8 Xl-^X 



T,._ 



exp \ 

Iminy 



0.65arcsec/ 



) , 



(A4) 



and 



(M,2p,±>(^)^4xl0"^x 
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Texp_\ f ^atm \ ^ ,^g^ 

Iminy \0.65arcscc/ 

Although the above relationships are quite crude, those 
are useful to quickly evaluate a magnitude of atmospheric 
PSF ellipticities for a given set of observational parame- 
ters. Note that these relationships are valid only for a tele- 
scope with a diameter of 8.2m. It may be worth pointing 
out that it is found from numerical simulations with differ- 
ent telescope diameters that the amplitude of PSF ellip- 
ticities is insensitive to D. This differs from a scaling rela- 
tionship of Pe oc D^ expected from the above consideration 
(likewise for {M^ j_)). This difference may be explained 
by the diameter dependence of the PSF which is not taken 
into account in the above consideration. According to 
the diffraction theory (e.g., Hardy 1998; Saha 2007), the 
PSF is the inverse Fourier transform of the optical transfer 
function (OTF) which is the auto correlation function of 
the product of the telescope pupil function with the atmo- 
spheric phase function. In that relationship, the diameter 
of telescope enters through the pupil function. Therefore, 
roughly and intuitively speaking, the power spectrum of 
PSF ellipticities relates to the OTF, and the diameter acts 
as a large-scale cutoff. Since the turbulence power spec- 
trum has a larger power on a larger scales (Kolmogorov 
power spectrum has the shape of /c"^^/"^), the large-scale 
cutoff imposed by the telescope diameter may have a great 
impact on the amplitude of PSF ellipticities. 



